Electronic and transport properties of azobenzene monolayer junctions as molecular switches 
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We investigate from first-principles the change in transport properties of a two-dimensional azobenzene 
monolayer sandwiched between two Au electrodes that undergoes molecular switching. We focus on trans- 
port differences between a chemisorbed and physisorbed top monolayer-electrode contact. The conductance of 
the monolayer junction with a chemisorbed top contact is higher in trans configuration, in agreement with the 
previous theoretical predictions of one-dimensional single molecule junctions. However, with a physisorbed 
top contact, the "ON" state with larger conductance is associated with the cis configuration due to a reduced 
effective tunneling pathway by switching from trans to cis, which successfully explains recently experimental 
measurements of azobenzene monolayer junctions. A simple model is developed to explain electron transmis- 
sion across subsystems in the molecular junction. We also discuss the effects of monolayer packing density, 
molecule tilt angle, and contact geometry on the calculated transmission functions. In particular, we find that a 
tip-like contact with chemisorption significantly affects the electric current through the cis monolayer, leading 
to highly asymmetric current-voltage characteristics as well as large negative differential resistance behavior. 

PACS numbers: 85.65.+h,31.15.E-,73.50.-h 
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I. INTRODUCTION 

Molecular junctions that incorporate photochromic 
molecules as reversible photo-switches between two differ- 
ent conductance states ("ON" and "OFF") have advanced 
considerably in the last decade-^. Azobenzene and its 
derivatives are the most frequently studied candidates for 
photoresponsive molecular switch, based on their confor- 
mational changes from a more thermodynamically stable 
trans configuration to a cis configuration in response to an 
external stimulus such as UV light^, and vice versa upon 
exposure to visible light or thermal excitation 4 . Previous 
first-principles studies for azobenzene molecular junctions 
focus on electron transport through ideal single-molecule 
junctions with one-dimensional (ID) electrodes and predict 
that the junction with trans configuration has a conductance 
higher than that with the cis configuration^—. 

More recent attention in experiments has been paid to the 
self-assembled monolayer (SAM) systems, which are highly 
ordered arrays of molecules on a two-dimensional (2D) sur- 
face with a chemisorbed bottom contact. A variety of experi- 
mental methods are used to apply a second contact on top of 
the azobenzene SAM and measure the photoinduced changes 
in its conductance, with sufficient flexibility to adapt to the 
height change of the SAM after its photo-switch, including 
liquid metal contact (eg., Hg drop£) and conducting atomic 
force microscopy'*'-'. These pioneering experiments found 
that the "ON" state with larger measured current is associ- 
ated with the cis configuration, contrasting previous theoreti- 
cal predictions based on single-molecule junction models. 

These experiments using 2D monolayer junctions present 
more complicated situations in which the ideal ID single- 
molecule junction model may not apply. First, in the mono- 
layer junction the two monolayer-electrode contacts are not 
symmetric. The one at the bottom interface is a chemisorbed 
contact, in which the end group of the molecule is chemically 
bonded to the electrode. The top contact is merely physical, 
i.e., no chemical bonding formed between the SAM and the 
top electrode. It is known that the way the molecule inter- 



acts with the electrode plays a crucial role in electronic trans- 
port through the molecule. The conductance of a junction 
with and without chemicontacts can differ by several orders of 
magnitud e) 11,12 . In addition, SAMs can form with differences 
in packing density and tilt angle, and in densely packed SAMs 
the conductance may involve both intramolecular and inter- 
molecular contribution due to molecules in close proximity 13 . 
Other effects such as the surface topography and roughness 
of the contact also influence the final geometry of the contact 
and thus lead to more complex transport properties. These un- 
certainties in the microscopic details of a monolayer junction 
make understanding the effects of monolayer-electrode cou- 
pling, intermolecular interaction and contact geometries crit- 
ical to elucidating its transport mechanisms. First-principles 
studies are especially valuable in this regard because informa- 
tion concerning the microscopic details of a monolayer junc- 
tion and their effects on electron transport are often impossible 
to obtain directly in experiments. 



Here we report a first-principles study of the electronic 
structure and transport properties of azobenzene monolayer 
junctions. The monolayer junction is constructed by attach- 
ing a Au( 111) electrode on top of an azobenzene monolayer 
chemisorbed on Au(l 1 1) substrate. We focus on clarifying the 
role of monolayer-electrode bonding in relation to the trans- 
port properties by investigating the junctions with two types 
of top monolayer-electrode contact: 1) a strongly bonded 
contact in which the molecules are covalently bonded to the 
Au surface, and 2) a weakly bonded contact which is ph- 
ysisorbed through van der Waals (vdW) interation. We find 
that both the trans to cis transformation of the molecules 
and the monolayer-electrode bonding play critical roles in de- 
termining the electron transmission as well as the current- 
voltage characteristics of the monolayer junction. The impact 
of monolayer packing density and tilt angle are discussed. We 
also consider the effect of contact geometry by adding an addi- 
tional Au atom at the top contact and find substantive changes 
in transport depending on this particular contact geometry. 



II. COMPUTATIONAL METHODS 

We use two density functional theory (DFT) based compu- 
tational methods in our calculations. First, all structural op- 
timizations, including calculations for both the azobenzene 
monolayer on a Au(lll) surface and the molecular junc- 
tion consisting of two monolayer-electrode contacts, are per- 
formed using the plane-wave-basis-set Vienna ab initio sim- 
ulation package (VASP) 14 . Then, electron transport calcu- 
lations of the molecular junction are performed using the 
atomic-orbital-basis-set TRANSIESTA code 15 , which incor- 
porates nonequilibrium Green's function formalism^ and 
DFT treatment of the electronic structure as implemented in 
the SIESTA package^. 

In the structural optimization calculations, projector aug- 
mented wave potentials with kinetic energy cutoff of 400 eV 
are employed. For the exchange and correlation functional 
we utilize the Perdew-Burke-Ernzerh generalized gradient ap- 
proximation (PBE-GGA)i£. A 3 x 3 x 1 fc-point sampling is 
applied based on the Monkhorst-Pack l9 scheme. All geome- 
tries are optimized until the remaining force on each ion fall 
below the convergence criterion of 0.02 eV/A. 

In the electron transport calculations, we utilize norm- 
conserving nonlocal Troullier-Martins pseudopotentials and 
the PBE-GGA exchange correlation functional, with a single- 
( plus polarization basis for the Au atoms and a higher level 
double-f basis set for the molecules. The equivalent plane- 
wave cutoff for the real-space grid of 150 Ry is used through- 
out the calculations. A 5 x 5 x 1 Monkhorst-Pack £-point 
sampling is used for self-consistent calculation, and a 9 x 9 
fc-point sampling is used for transmission calculations. The 
transmission function is obtained by averaging transmission 
coefficients in each fc-point for an applied bias voltage V, 
T(E, V) = £iUL T klA (E, V), where 



7U, = Tr[r 1 G R r 2 G A ] 



(1) 



is the ^-resolved transmission coefficient, G R ' A are the retarded 
(advanced) Green's function matrices, and T^ = «[£i,2 - £{ 2 ] 
where the Z12 are the retarded self-energies due to the exis- 
tence of bottom (top) semi-infinite electrodes. The electric 
current as a function of the applied voltage is obtained by the 
integration of the transmission function 



I(V) 



- 2 il 



T(E,V)[f(E-m)-f(E-^ 2 )]dE. (2) 



where e is the elementary charge, / is the Fermi-Dirac distri- 
bution function, pt\2 — Ep ± eV/2 are the chemical potentials 
for the bottom (top) electrodes. 



III. AZOBENZENE MOLECULAR MONOLAYER ON 
AU(lll) SURFACE 

A. Structural models 

For the azobenzene monolayer on the Au(lll) substrate 
(hereafter we denote this as Au-AB), the systems that we 
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FIG. 1 : (Color online). Equilibrium configurations of trans (a) and 
cis (b) azobenzene molecules chemisorbed on Au(l 1 1) surface via a 
CITS linker. Insets: side view of the corresponding geometry. C, 
H, N, S and Au atoms are colored gray, green, red, blue and yellow, 
respectively. 



study consist of an azobenzene monolayer chemisorbed on 
the surface of Au(lll) substrate by a CH2S linker for each 
molecule, as illustrated in Fig.Q~]produced by VESTA graphi- 
cal program 20 . The substrate is modeled by a slab consisting 
of five Au(l 1 1) layers. The molecular coverage is kept fixed 
at one molecule per ( V3 x y/3)R30° surface unit cell, which 
corresponds to a packing area of 90.8 A per molecule. The 
supercell thus consists of 87 atoms comprising 12 atoms per 
Au layer and 27 atoms for the azobenzene molecule with the 
linker. A large vacuum spacing of 15 A is used between the 
topmost molecular atom and the next Au slab to prevent in- 
teraction between adjacent images. All atoms except the two 
bottom Au layers are fully relaxed. 



B. Adsorption geometry and energetics 

In order to determine the preferred adsorption position for 
the azobenzene monolayer on the Au(l 1 1) surface, three pos- 
sible initial binding sites of both azobenzene isomers with the 
bottom CFFS linker are optimized on the Au(l 1 1) surface, the 
on-top site, the fee- and hep-bridge sites, as shown in Fig. [2] 
The calculated total energy and binding energy of an azoben- 
zene molecule at different adsorption sites are reported in Ta- 
ble U For both trans and cis configurations, we find the most 
stable adsorption position to be the fee -bridge site. We also 
find the trans configuration to be more stable by about 0.6 
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FIG. 2: (Color online), (a) Top view indicating the lateral position 
of the sulfur atom in the structural optimizations carried out to de- 
termine the most stable adsorption site of the azobenzene molecule 
with CH 2 S linker on Au(lll) surface, (a): initial position; (b): the 
optimized position. 



TABLE I: Calculated total energy (in eV per molecule), binding en- 
ergy (in eV per molecule) and charge transfer (in |e|) of an azoben- 
zene molecule on Au(lll) surface at different adsoiption sites. The 
total-energy of Au/tram AB at fcc-bridge site is set to be reference 
zero. 
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FIG. 3: (Color online). Total energy as a function of tilt angle for 
trans azobenzene molecules chemisorbed on Au(lll) surface. The 
left and right insets are side views of the corresponding optimized 
structures with a tilt angle of 20° and 45°, respectively. 



monolayer is usually determined by other factors, such as am- 
bient temperature and monolayer packing density. It has been 
found that the molecular monolayer undergoes a phase tran- 
sition from a tilted structure to a vertical structure at room 
temperature 24 . Also, in densely packed monolayers at full 
coverage, monolayers with a small tilt angle can be generally 
expected, with molecules standing with their long molecular 
axes close to the surface normal. 



eV than the cis configuration, in good agreement with the re- 
sults from previous calculations 21 and experiments 2 .. Bader 
analysis based on the real-space charge density^ for the Au- 
AB shows a small amount of charge transfer from Au sub- 
strate to the molecule for all cases, as given in Table U The 
charge transferred to the molecule is primary localized at the 
CH2S linker, and not efficiently shared with the core part of 
the molecule. Hereafter we present only the results for sys- 
tems with the azobenzene monolayer chemisorbed on the fcc- 
bridge site of the Au(l 1 1) surface. 

We also optimize the tilt angle of the azobenzene mono- 
layer on Au(lll) surface. In Fig. [3] we show the calculated 
total energy of a trans azobenzene monolayer on Au( 1 1 1 ) as a 
function of tilt angle ranging from 10° to 45°. At every single 
energy point in Fig. [3] the monolayer structure is optimized 
by performing geometry relaxation with an initially given tilt 
angle. Within this angle range, we find that the lowest en- 
ergy structure is the 20° tilted molecular chemisorbed on the 
Au(l 1 1) surface. We expect the same finding for the cis con- 
figuration since the trans to cis transformation only change 
the geometry of top half of the azobenzene molecule. There- 
fore, for the rest of the study we choose azobenzene mono- 
layer with tilt angle of 20° as the equilibrium structure for 
both trans and cis configurations, as shown in Fig. Q] 

However, it should be noted that in Fig. [3] the calculated 
total energies of these systems with different tilt angles are 
very close. In experiments, the tilt angle of molecules in the 



C. Electronic structure 

Figure |4] analyzes the interaction between the azobenzene 
and Au( 111) surface in terms of the density of states projected 
onto the atoms (PDOS) in the molecule as well as the sul- 
fur linker. For both trans and cis configurations, we see the 
two frontier orbital levels of azobenzene, highest occupied 
molecular orbital (HOMO) and lowest unoccupied molecu- 
lar orbital (LUMO), lie at each side of the Au Fermi level. 
Unlike the other molecular orbital levels nearby, the HOMO 
and LUMO contain large contributions from the nitrogen pair 
of the azobenzene. The position of HOMO in the cis config- 
uration is closer to the Au Fermi level than that in the trans 
configuration. This is quite different compared to the results 
from a study of a single azobenzene molecule on an Au(l 1 1) 
surface 21 , in which the downshift of molecular orbital level is 
more significant for the cis isomer, resulting in the HOMO of 
cis azobenzene being away from the Au Fermi energy and its 
LUMO being very close to the Fermi energy. This difference 
arises because the azobenzene molecules in the monolayer are 
standing upright (nearly perpendicular to the surface), result- 
ing in a weak interaction between the molecule and the bottom 
surface, whereas a single azobenzene molecule orients paral- 
lel to the surface and thus the closer proximity between the 
azo group (N pair) of the molecule and the Au surface creates 
a stronger interaction between the m 21 ' 25 . For trans monolayer 
with a tilt angle of 45° on the Au surface, we find essentially 
the same features in the PDOS as obtained for 20° tilted trans 
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FIG. 4: (Color online). Projected density of states of the molecule 
and various atoms of trans (a) and cis (b) azobenzene monolayer 
chemisorbed on Au(lll) surface. For trans configuration, addition- 
ally shown is the PDOS of the molecule with a tilt angle of 45° on 
the Au surface. Vertical green line represent the Fermi energy which 
is at eV. 



monolayer, but all molecular orbital levels are roughly equally 
shifted by about 0.4 eV to lower energies. 



IV. AZOBENZENE MONOLAYER JUNCTIONS AND 
TRANSPORT PROPERTIES 

A. Structural models and electronic structures 

Once the equilibrium geometry of the Au-AB is obtained, 
we construct the corresponding molecular junction by extend- 
ing the bottom Au(l 1 1) substrate into a six atomic layer slab, 
and attaching a second Au(lll) slab consisting another six 
atomic Au layers on top of the azobenzene monolayer. To ex- 
amine the effect of monolayer-electrode contact, two types of 
interfacial contacts for the top electrode are studied: (1) the 
top Au surface is covalently bonded to the azobenzen mono- 
layer with CH2S linkers same as those at the bottom contact 
(hereafter we denote this type of junction as Au-AB-Au), and 
(2) the top contact is physisorbed on the azobenzene mono- 
layer surface without the CH2S linkers (we denote this as Au- 
AB Au junction). We optimize the top interfacial contact by 
computing the total energies of the system as a function of 
the distance between the left and the right electrodes. Every 
single energy point is calculated by performing geometry op- 
timization with a constrained electrode-electrode separation. 
Therefore the equilibrium geometry is obtained as the dis- 
tance at which the total energy is minimal. Due to the failure 
of DFT GGA exchange correlation functionals to account for 
vdW interactions, we optimize the electrode-electrode sepa- 
ration using vdW density functionals optB88-vdW— for junc- 
tions with physisorbed contacts. Fig.|6](a) and (b) shows the 
optimized monolayer junction geometries with equilibrium 
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FIG. 5: (Color online). Projected density of states of the molecule in 
trans (a) and cis (b) azobenzene monolayer junctions. 



electrode-electrode distance for trans and cis configurations, 
respectively. 

We show the PDOS of the azobenzene molecules in the 
monolayer junctions in Fig. [5] As compared to the PDOS 
of the bare monolayer on Au(lll) surface, it can be clearly 
seen that all molecular orbital levels, including the HOMO 
and the LUMO, are shifted to lower energies with the pres- 
ence of a top monolayer-electrode contact. In Fig. |5](a), the 
downshift of molecule orbital levels in the Au-?ranjAB|Au 
junction is only minimal, as compared to Fig.|4](a), primar- 
ily because of a weak bonding between the top electrode and 
the transAB molecule. In contrast, with the sulfur linkers 
in the top contact of the Au-fransAB-Au junction, a stronger 
monolayer-electrode bonding significantly lower all molecu- 
lar orbital levels in energies by about 0.3 eV. The downshift 
is even more significant in the cases of cis configurations as 
shown in Fig.[5](b). All cis azobenzene orbital levels shift to 
lower energies by about 0.4 eV in Au-c is AB Au and 0.6 eV in 
Au-c/jAB-Au, compared to Fig.|4](b). 



B. Transmission functions at zero-bias 

In Fig. 0(c) and (d) we show the calculated transmission 
functions at zero-bias for azobenzene monolayer junctions 
with trans and cis configurations, respectively. For the Au- 
transAB-Au junction with chemisorbed top contact, transmis- 
sion function has three broad resonance peaks with near per- 
fect transmission amplitudes shown in the given range of en- 
ergy from -2 eV to 2 eV related to the Fermi energy of the 
Au electrodes: two are below the Fermi energy (about - 1 .6 eV 
and -0.7 eV), and the third is above it (about 1.1 eV). These 
three transmission peaks can be directly associated to the fron- 
tier molecular orbitals HOMO-1, HOMO and LUMO, respec- 
tively, despite some minor changes in the energy as compared 
to the PDOS in Fig. [5] due to the periodic boundary condi- 
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FIG. 6: (Color online). Equilibrium geometries (a-b) and zero-bias electron transmission functions (c-d) of trans and cis monolayer junctions 
with two types of top contact. In panels (c) and (d) the y-axis is in logarithm scale, and the Fermi level is shifted to 0. Ad represent the change 
in the electrode-electrode separation for a junction with its geometry slightly away from the equilibrium structure. 



tion used in the supercell calculation for PDOS. The trans- 
mission peak associated with the HOMO of trans azobenzene 
is relatively narrower in width and lower in height because 
of the fact that the HOMO is more localized in the N pair of 
the molecule and thus contributes less to the electron trans- 
mission. These resonance peaks in the transmission function 
are away from the Au Fermi level, and zero-bias conductance 
of the monolayer junction is primarily dominated by non- 
resonance electron tunneling at the Fermi level, with a trans- 
mission coefficient T ~ 0.02. For the Au-fra«iAB| Au junc- 
tion with a physisorbed top contact, we find that the calculated 
transmission function is significantly decreased by about two 
orders of magnitudes as compared to that of Au-fra/wAB-Au 
junction, and, at the Fermi level, T ~ 0.0004. The resonance 
transmission peaks of Au-fra«jAB|Au junction shift about 0.2 
eV upward in the energy, which follows closely the peaks in 
PDOS as shown in Fig.|4](a). Besides three major peaks, there 
are several small resonance peaks appearing in the transmis- 
sion function of Au-fra«iAB|Au around eV to 1.2 eV which 
are induced by the vacuum gap in the top contact. Their rel- 
atively smaller amplitudes provides only finite contribution to 
the tunneling current under certain bias voltages, which we 
will discuss later in the next section. 

For monolayer junctions with cis configuration, however, 
removing the CH2S linkers in the top contact leads to a rel- 



atively much smaller decrease in transmission function near 
the Fermi energy, as shown in Fig. |6](d). The transmission 
functions at the Fermi level for junctions Au-cwAB-Au and 
Au-cwAB|Au are T w 0.006 and 0.002, respectively. For 
both Au-c/sAB-Au and Au-c/jAB|Au, within the given range 
of energy two resonance peaks can be found in the transmis- 
sion functions, which are associated with the cis azobenzene 
HOMO and LUMO. They are also away from the Au Fermi 
level, and the peaks of Au-cwAB-Au junction are lower in en- 
ergy than that of Au-c«AB|Au, corresponding to the PDOS 
peaks as shown in Fig.|4](b). 

The above calculated zero-bias transmission function can- 
not be simply studied by its relation to the electrode-electrode 
separation or the azobenzene molecular length. Instead, we 
can describe the total transmission function that reflects the ef- 
ficiency of electronic transport from one electrode to the other 
through the azobenzene monolayer junction as 



" c mol mol tc 



(3) 



where T\, c and T tc give the efficiency of electron transport 
across the bottom and top contacts, and T p , and T p , re- 

r ' mol mol 

fleet the electron transport through the first (bottom half) and 
second (top half) parts of the azobenzene molecule, respec- 
tively. From Fig. [6] (a) and (b), it is safe to assume that the 
first term T\, c in the above transmission equation is the same 



for all four junctions. However, the other three terms, T p v 

7^ ol and 7^, are quite different from case to case, result- 
ing in different total transmission coefficients at the Fermi 
energy for the above junctions. In the case of the sym- 
metric Au-transAB-Au junction, we have 7 P , — 7 P , and 

J ' mol mol 

The - T tc , and the junction has a large transmission coef- 
ficient at the Fermi level, 7 ~ 0.02. The absence of the 
CH2S linker in the Au-transAB\Au junction causes a large 
decrease in T tc , and reduces the total transmission coefficient 
to 7 ~ 0.0004. One may expect such a large change in the 
transmission term 7 tc to also apply to the cis configuration 
junction. However, as shown in Fig.|6](d), the total transmis- 
sion coefficient for Au-cwAB|Au (7 ~ 0.002) is only about 
three times smaller than that for Au-cisAB-Au (T ~ 0.006). 
This clearly indicate that, in the Au-c/*AB|Au junction, the 
transmitted electron travel directly from the bottom contact to 
the top contact through only the first (bottom) part of the cis 
azobenzene molecule, and avoid the second (top) part of the 
molecule. Thus for the Au-cwAB|Au junction the effective 
tunneling pathway is greatly reduced, resulting a total trans- 
mission function T - 7 DC ■ 7^ ]ol ■ 7 tc . Note this is equal to Eq. 

|3]with term 7^ = 1 that compensates the smaller term 7 tc of 
the Au-c«AB|Au junction, resulting in a smaller difference in 
the total transmission coefficients between the Au-cwAB|Au 
and Au-cisAB-Au junctions as compared to that of the junc- 
tions in the trans configuration. 

It should be also noted that the transmission terms 7 P , and 

mol 

7 P , are also different between junctions with trans and cis 

mol J 

configurations. Non-resonant electron transmission through 
a molecule 7 mo i can be empirically expressed as 7 mo i = 
exp(—pt), where / is the length of the molecule and represents 
the width of effective tunneling barrier, and /? is the tunneling 
decay factor given by /J = (2 ~J2m*(pa)/n where cf> is the barrier 
height for tunneling that is determined by the frontier molec- 
ular orbital level related to the Fermi energy, <p = Ep - E M o, 
m* is the effective mass and a is the shape constant of the bar- 
rier. For junctions with chemisorbed top contact, at the Fermi 
energy the total transmission ratio between trans and cis con- 
figurations is t = T'™ s /T™ } = rapRPcfa -PtransWl Giving 
t ~ 3.3 from Fig. [6] and / m 4.5A estimated from the length 
of each half of the azobenzene molecule (remain unchanged 
upon switching), we have/?,^ — fitrcms - 0.13A -1 . For junc- 
tions with physisorbed top contact, from the above discussions 
we know that in the Aa-trans AB | Au the effective barrier width 
for the expression of T'^" s is 21 and that in Au-cisAB|Au is /, 
thus the total transmission ratio between these two junctions 
at the Fermi energy should be t' = exp[(J3 C j S l - f3 trans 2l)]. Giv- 
ing t' x 0.2 from Fig. [6] we get the tunneling decay factors 
p trans = 0.49 A -1 andjScjj = 0.62A"1. Both of them are in 
good agreement with the literature values (J5 m 0.4 - 0.6A -1 ) 
commonly obtained for short ^-conjugated molecules 2 ^. A 
possible cause for a smaller tunneling decay factor f3, rans than 
fids could be the fact that the HOMO level of junction in 
the trans configuration lie closer to the Fermi level than the 
HOMO level of the junction in the cis configuration, as shown 
in Fig. [5] and Fig.|6j resulting in a lower tunneling barrier, <p. 
The barrier shape constant a could also be affected after the 



switching of azobenzene monolayer from trans to cis form. 

In experiments, the monolayer-electrode contact of the 
junction may not be exactly at its equilibrium position, but 
rather be slightly stretched or compressed. To simulate this 
situation, we calculate the transport properties of junctions 
with an increased (or decreased) electrode-electrode separa- 
tion Ad = 0.5A (or Ad = -0.5A) with regard to its equilib- 
rium value (Ad = 0) for both trans and cis monolayer junc- 
tions. The transmission results for stretched and compressed 
junctions are also shown in Fig. [6] We find that the conduc- 
tance of the Au-AB-Au junction with chemisobed top contact 
is insensitive to the change of Ad for both trans and cis config- 
urations, where as for Au-AB|Au junctions with physisobed 
top contact, a small change in the electrode-electrode separa- 
tion in either direction changes the conductance significantly. 
This demonstrates that, with a physisorbed top contact, the to- 
tal transmission function of the monolayer junction is highly 
related to the 7 tc (contribution from the contact) which is very 
sensitive to the monolayer-electrode contact distance. In the 
case of Au-cisAB|Au junctions, we also notice that the two 
peaks in the transmission function at about -0.7 eV and 1.1 eV, 
which are associated with the HOMO and the LUMO, move 
up in energy with increased electrode-electrode separation. 
This phenomenon has also been noticed for other stretched 
molecular junctions 28 , in which the shift of HOMO resonance 
in the transmission function toward the Fermi level increases 
the transmission coefficient at the Fermi energy. However, 
here the decrease in the transmission coefficient at the Fermi 
level, due to the stretching of the junction, is much more 
significant. The transmission peak associated with the cis 
HOMO is still more than 0.5 eV away from the Fermi level 
with Ad = 0.5A, making the decrease in low bias conduc- 
tance unlikely to be compensated. In contrast, for stretched 
and compressed Au-fran*AB|Au junctions we do not see such 
obvious change in the positions of transmission resonances. 

In order to understand the possible influences of the mono- 
layer tilt angle 9 on the transport properties, we carry out fur- 
ther calculations for Au-trans AB-Au junction with different 9. 
The optimized junction structures with 45° or 10° tilted trans 
monolayer are shown in Fig. |7] (a) and (b), respectively. A 
large change in the electrode-electrode separation occurs be- 
tween these two cases (Ad = 2.5A). However, as shown in 
Fig- 0(c), for junctions with a 45° or 10° tilted trans mono- 
layer, we find essentially the same features in the transmission 
functions as those obtained for the optimized geometry with 
9 = 20°, yet with a small change in position of transmission 
resonance peaks associated with the trans azobenzene molec- 
ular orbital levels for 9 = 45°. The transmission at the Fermi 
level for 45° tilted junction is only slightly increased as com- 
pared to the junctions with small tilt angles. The source of 
this increase is the downward shifting of the broad peak in the 
transmission function associated with the trans LUMO. The 
shifting of molecular frontier levels resulting from the change 
in tilt angle may be resulted from a change in the effective 
dipole moment as well as the work-function at the molecule- 
metal interfaces 2 ^. Similar phenomena have also been noticed 
in tilted alkanethiol monolayers^, in which the tilt brings the 
HOMO resonance closer to the Fermi energy. Nevertheless, 
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FIG. 8: (Color online). Current voltage characteristics of Au-AB-Au 
(a) and Au-AB|Au junctions with trans and cis configurations. 



FIG. 7: (Color online). Equilibrium geometries (a-b) and zero-bias 
transmission functions for Au-frarcsAB-Au junctions with a tilt angle 
(9of45 o andl0 o . 



we conclude that for dilute monolayer junctions without in- 
termolecular contribution, the molecular tilt has very little ef- 
fect on its conductance since the effective tunneling pathway 
through the molecule remains unchanged upon tilting. 



C. Current-voltage characteristics 

The self-consistently calculated current-voltage character- 
istics (I-V) of Au-AB-Au and Au-AB|Au junctions are de- 
picted in Fig.|8](a) and (b), respectively. For the Au-fransAB- 
Au junction, a highly symmetric I-V curve arises as expected 
with two symmetric monolayer-electrode contacts, and its cur- 
rent is larger than that of the Au-c/sAB-Au junction in small 
bias region, consistent with the zero-bias transmission func- 
tions. Thus, with chemisorbed top contact, the two isomers 
realize different conductance states of the monolayer junction 
corresponding to "ON" {trans configuration) and "OFF" (cis 
configuration). The current of the Au-fransAB-Au junction 
becomes smaller than that of the Au-c/sAB-Au only in the re- 
gion of large negative bias voltage (-1.2 V to -2.0 V). This 
can be attributed to the fact that the Au-c/sAB-Au junction 
is asymmetric in geometry, resulting in an asymmetric I-V 
curve which has a slightly larger current in the large negative 
bias region. Such asymmetric behavior is also shown in the 
I-V curves for junctions Au-fransAB|Au and Au-c/.sAB|Au, 
with asymmetric top (physisorbed) and bottom (chemisorbed) 
contacts. Nonetheless, for junctions with a physisorbed top 
contact, the cis configuration ("ON" state) exhibits a larger 
current for any given voltage than that of the trans configu- 



ration ("OFF" state), and the ON/OFF ratio of the current is 
more pronounced than that of the junctions with chemisorbed 
top contacts. These theoretical results for junctions with a ph- 
ysisorbed top contact are, therefore, in reasonable qualitative 
consistency with the reported measurements^^. A quantita- 
tive matching of the ON/OFF ratio to the experiments is not 
expected because the length of the molecules in experiments 
are much longer than those in our study. 

In Fig. [8] (b), we find the current of the Au-fran.yAB|Au 
junction oscillates as a function of positive bias voltage, which 
correspond to the small resonance peaks found in the trans- 
mission function above the Fermi level as shown in Fig.|6](c). 
As we increase the voltage, these resonance peaks sequen- 
tially join the bias window and contribute to the current, re- 
sulting in resonant tunneling current at certain bias voltages. 
However, these resonances can only provide finite contribu- 
tion to the current in the trans configuration and the current 
ratio between the trans and cis configurations is almost unaf- 
fected. 



D. Effect of molecular packing density 

The above calculations are based on a diluted monolayer 
structure with a low packing density, in which the intermolec- 
ular bonding between molecules in the monolayer can be con- 
sidered negligible and junction transport properties are pri- 
marily dominated by the intramolecular contribution. To in- 
quire into effect of intermolecular interaction on the transport, 
Au-AB-Au junctions with a higher monolayer packing den- 
sity are also considered, in which the molecules are positioned 
closer to each other to increase the molecular interactions. In 
this case for each junction we have two azobenzene molecules 
per ( V3 x V3)/?30° Au(l 1 1) unit cell, and the packing area 



8 



(a) 



Au-transAB-Au 




Au-2fra/isAB-Au 



*' ft' «■ 



(b) Au-c/'sAB-Au 




Au-2c/'sAB-Au 




Energy (eV) 

FIG. 9: (Color online). Equilibrium geometries (a-b) and zero-bias 
transmission functions (c) of azobenzene monolayer junctions with 
different packing density. The calculated transmission coefficients 
for junctions with high monolayer packing density are divided by two 
(number of azobenzene molecules in one cell) in order to compare 
with the results with low packing density. 



corresponds to 45.4 A per molecule. The optimized junction 
structures for trans and cis configurations are shown in Fig. [9] 
(a) and (b), and we denote them as Au-2fra«.yAB-Au and Au- 
2c/sAB-Au, respectively. For the cis configuration, as com- 
pared to the Au-cwAB-Au junction with a low packing den- 
sity, we see that the azobenzene molecules in the monolayer 
of the Au-2cwAB-Au junction become distorted, so that the 
top half of the molecule is no longer parallel to the top Au 
surface but instead rotated relative to the Au surface, resulting 
from an increased intermolecular interaction with decreased 
spacing between the molecules. 

The calculated zero-bias transmission functions are shown 
in Fig. |9](c). For the trans configuration, doubling the pack- 
ing density produces only negligible changes in transmission 
function. However, for the cis configuration, the junction's 
transmission coefficients are significantly reduced especially 



in the energy region around the Fermi level. The decrease in 
conductance by increasing the packing density has also been 
observed in other monolayer junctions 13 , but the reason has 
been attributed to a shift of the LUMO resonance peak in 
the transmission function to a higher energy away from the 
Fermi level. Our result differ from the Ref. [lj in that the 
LUMO peak in the transmission function of the Au-2c;\sAB- 
Au junction slightly shifts to a lower energy which is closer to 
the Fermi level as compared to the Au-c/sAB-Au. Neverthe- 
less, the transmission coefficients for the HOMO resonance 
peak and for the energy region around the Fermi level are de- 
creased by about one order of magnitude. The considerable 
decrease in the transmission function can be attributed to the 
geometry distortion in the cis azobenzene monolayer which is 
caused by intermolecular interactions with the dense packing, 
as shown in Fig.|9](b). Our results show clearly that interac- 
tion among the azobenzene molecules of the monolayer junc- 
tion can induce large changes in transport in the cis configu- 
ration, hence suggesting that higher conductance ON/OFF ra- 
tio can be achieved in densely packed azobenzene monolayer 
junctions with chemisorbed monolayer-electrode contact. For 
junctions with a physisorbed contact, since the "ON" state is 
associated with the cis configuration, we expect that increas- 
ing the monolayer packing density could decrease the conduc- 
tance ON/OFF ratio. 



E. Effect of top-contact geometry 

Having addressed the transport properties of azobenzene 
molecular junctions with a perfect contact surface at the 
monolayer-electrode interface, it is necessary to also study the 
consequences of non-ideal situations. To investigate the role 
of contact atomic structure, we consider a simple but possible 
situation: the sulfur atom bound to a single Au atom protrud- 
ing from the surface of the top Au electrode, which is shown 
in Fig.[10](a) and (b) for trans and cis configurations, and we 
denote them as denoted as Au-AB-Aui-Au and Au-2c;\sAB- 
Au, respectively. The purpose of this consideration is to sim- 
ulate possible situations in experiments in which the surface 
of top electrode may not be atomically flat when contacting 
the molecules in the monolayer. 

The calculated transmission functions are shown in Fig.flOl 
(c). We find that the additional Au atom decreases the trans- 
mission coefficients at non-resonance tunneling regions in the 
transmission function for both trans and cis configurations, 
as compared to Fig. [6] Moreover, the resonance peaks in the 
transmission functions associated with the frontier molecular 
orbitals, all shift up in energy at a noticeable change. The 
HOMO transmission resonance of Au-fran^AB-Aui-Au junc- 
tion shifts from -0.7 eV to -0.4 eV; while that of Au-dsAB-Au 
junction shifts from -0.85 eV to -0.2 eV, which is very close 
to the Fermi energy. 

In Fig. [10] (d) we show the calculated current- voltage char- 
acteristics for Au-AB-Aui-Au junctions. The most noticeable 
feature appearing in Fig. [10] (d) is the asymmetry in the cis 
configuration. For a negative bias voltage below 1 eV, the cur- 
rent through the cis configuration is significantly higher than 
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FIG. 11: (Color online). Electron transmission coefficients of Au- 
AB-Aui-Au junctions with trans and cis configurations under vari- 
ous bias voltages from -2.0 to 2.0 V. The region between two green 
vertical lines indicates the bias window. The arrow indicate the bias- 
induced transmission peak in the transmission function in the cis con- 
figuration. 



FIG. 10: (Color online). Equilibrium geometries (a-b), zero-bias 
electron transmission coefficients (c) and current-voltage character- 
istics of Au-AB-Au r Au junctions. 



the current through the trans configuration. As the negative 
bias is further increased, the current through the cis configu- 
ration decreases dramatically with increasing voltage, result- 
ing in a region of negative differential resistance (NDR). This 
NDR feature also appears in the positive bias region (from 
1.2 V to 1.7 V) in the cis configuration, but is not present for 
any bias region in the trans configuration. In contrast to Au- 
c/sAB-Aui-Au, the I-V characteristics of Au-transAB-Au\- 
Au junction are basically unchanged by the presence of the 
additional Au atom as compared to Fig.[8](a). 

To shed light on the origin of the NDR effect, we give the 
transmission function T(E, V) at a series of biases from -2.0 
to 2.0 V at 0.5 V intervals as shown in Fig. QT| In the case 
of the cis configuration, a voltage of -0.5 V creates a bias- 
induced transmission peak at about -0.2 eV (marked with 
an arrow in Fig. Q~l} adjacent to the resonance peak associ- 
ated with molecular LUMO. As the negative bias is further 
increased, this transmission peak, with a broad width, moves 



up in energy and towards the molecular LUMO resonance. It 
dominates the transmission function within the bias window, 
thus raising a much larger current than that at positive bias 
and resulting in a noticeable rectification effect for Au-c/sAB- 
Aui-Au junction in Fig. [I0](c). In the case of trans configu- 
ration, however, only several bias-induced transmission peaks 
exist at positive bias, and all of them have much lower am- 
plitudes and shorter width leading to only small current con- 
tributions. Correspondingly, in Fig. [10] (c) the rectification 
behavior in the I-V curve of Au-fra/wAB-Aui-Au junction is 
not obvious and there is no NDR effect. This demonstrate that 
the azobenezene monolayer junction can be used for photo- 
controlled molecular rectifier with proper engineering of its 
contact geometry. 

It should also be noted that for both trans and cis configu- 
rations, the bias drives the LUMO resonance peak away from 
the Fermi level, so its contribution to the current becomes 
important only at a sufficiently large bias (e.g. V ~ 2 V) 
which should be comparable to the HOMO-LUMO gap of the 
molecule. This is also true for the calculated I-V characteris- 
tics in Sec. |IV Cl Therefore, at low bias the electronic trans- 
port of all the monolayer junctions in consideration is mainly 
facilitated by non-resonant electron tunneling. 
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V. CONCLUSION 

Our first-principles calculations indicate that monolayer- 
electrode bonding, intermolecular interaction, and contact 
geometry all play very important roles in determining the 
conductance ratio of azobenene monolayer junction between 
trans and cis configurations. With a strongly bonded top con- 
tact from chemisorption, the zero-bias transmission of trans 
monolayer is higher than that of the cis monolayer. Replac- 
ing it by a weakly bonded contact with physisorption will de- 
crease the conductance by about two orders of magnitudes for 
the trans monolayer. However, in the case of the cis config- 
uration, the physisorbed contact shortens the effective tunnel- 
ing pathway, leading to a conductance that is less sensitive to 
weak bonding. The mechanism are interpreted by applying a 
simple transmission model to the calculated zero-bias trans- 
mission functions, which accounts for electron transmission 
through each subunit of the molecular junction and provides a 
clear physical picture for understanding the junction. The cal- 
culated current- voltage characteristics indicate that, for junc- 
tions with physisorbed top contact, the "ON" state with larger 
current is associated with the cis configuration, which are in 



agreement with recent experiments. Our calculations demon- 
strate the intermolecular interaction in the trans configura- 
tion is negligible even in a densely packed monolayer; while 
for cis configuration increasing the monolayer density causes 
a considerable distortion in the monolayer structure, result- 
ing a decrease in its transmission function thus increasing the 
conductance ratio. We also find that the calculated current- 
voltage characteristics for the trans monolayer remain sym- 
metric upon changing of the top contact geometry. In con- 
trast, a slightly modified contact geometry will affect the elec- 
tric current through the cis monolayer, leading to a highly un- 
symmetrical current-voltage curve as well as a large negative 
differential resistance behavior. These results suggest that the 
molecule-lead contact, the state of the molecule, the morphol- 
ogy of the metal surface, and the packing density of the mono- 
layer are all parameters in the play. Our investigations thus 
will deepen understanding electron transport through azoben- 
zene monolayer junctions. 
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